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Abstract 

We study the conditions for the existence of neutrino oscillations in the 
field-theoretical approach which combines neutrino production and detection 
processes in a single Feynman graph. The "oscillating neutrino" is represented 
by an inner line of this graph where, due to the macroscopic distance L 
between source and detector, the neutrino propagators for neutrinos with 
definite mass are replaced by the projection operators unto the neutrino states 
on mass shell. We use as a concrete model reaction the neutrino source and 
detector as given in the LSND experiment and we carefully take into account 
the finite lifetime of the stopped muons which provide the z?« beam. We 
show that the field-theoretical approach provides a solid method to locate 
all possible conditions and allows to separate unambiguously their different 
origins. Some of these conditions are independent of L whereas others state 
that coherence is lost when L exceeds a certain "coherence length". Also it 
turns out that, at least in the concrete situation considered here, the concept 
of neutrino wave packets is not supported by the field-theoretical approach for 
realistic experimental conditions, i.e., the neutrino energy spread is incoherent 
in origin. 

14.60.Pq, 03.65.-w 
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I. INTRODUCTION 



The standard treatment of neutrino oscillations provides a beautiful and simple pic- 
ture of this important phenomenon. With the mixing matrix relating the left-handed neu- 
trino flavor fields with the left-handed neutrino mass eigenfields defined by v La = £V U a jV L j 
(a = e, fi, t, . . .) it allows to derive the oscillation probabilities for antineutrinos 



P- 



x / m 2 L\ 

j v 

{( Am 2 L\ I 
U* pj U aj Up k U* ak exp ^-i—^-j I , (1.1) 



valid in the ultrarelativistic limit with Am^ = m 2 — m 2 k where ni\ < ni2 < . . . denote 
the neutrino masses, L is the distance between neutrino source and detector and E v is the 
neutrino energy. The probability for neutrinos is obtained from Eq. (|l.l| ) by the substitution 
U —>[/*. In the following, Greek indices always indicate neutrino flavors and Latin indices 
mass eigenstates or fields. However, after a closer look one discovers that the standard 
derivation of Eq. fll.lj ) needs clarification in several points (see, e.g., Ref. || for a summary 
of these problems). This has first been attempted by using neutrino wave packets @-12 



whereas Ref. has pioneered the idea of considering the complete neutrino production - 
detection chain using only those quantities for the description of neutrino oscillations which 
are really observed or manipulated in oscillation experiments |13|-[T5| in order to obtain 
unambiguous results. See also Ref. []16| for a sort of combined field theory - wave packet 



approach. The present interest in theoretical treatments of neutrino oscillations can be 
phrased by the following question: Under which conditions is formula (\1-J\) valid? Since the 
neutrino wave packet formalism does not work with the physical observables we find the field- 
theoretical approach treating neutrino production and detection |T3[ the most appropriate, 



unambiguous and general way to analyse the problem of coherence in neutrino oscillations. 
In particular, when several quantities defining a length are involved, an improvement of 
the wave packet approach is called for in order to distinguish the roles and origins of these 
lengths. 

The transition probability ( |1.1| ) is given by the square of the sum over the amplitudes of 
the neutrino mass eigenstates, i.e., by a coherent summation over the mass eigenstates. The 
first term in the second line of Eq. ( |L1| ) represents the purely incoherent summation over the 
mass eigenstates whereas the second term denotes the interference terms. The exponentials 
exp(—i27rL/L° s k c ) with the oscillation lengths defined by 



'ST - (i-2) 



show the oscillatory behaviour of the transition probability as a function of L/E v . Eq . ( |1 . 1|) 
is a theoretical expression without regard to an actual experimental situation. In the de- 
scription of a neutrino oscillation experiment it is possible that, after taking into account 
the experimental conditions, some or all of the interference terms drop out as a consequence 
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of certain averaging or suppression mechanisms to be discussed in the following. Note that 
the effect of such mechanisms is equivalent to a partial or complete incoherent summation 
over the neutrino mass eigenstates. 

One such class of mechanisms is given by all effects leading to an energy spread of the 
neutrino beam. It has been shown in Ref. J?] that if we label such effects by a then each of 
these effects giving an energy spread AE a leads to a coherence length 

t coh rose v / 1 q\ 

n a;jk ~ ^jk K 1 - 6 ) 

independent of the fact whether this spread has to be interpreted as a coherent or inco- 
herent effect. In the context of a neutrino energy spread "incoherent" means that single 
neutrinos have a definite energy but the neutrino beam has an energy spread whereas by 
"coherent" it is understood that a single neutrino state is a superposition of different en- 
ergies.Q Note, however, that in the field-theoretical approach the notions "incoherent" and 
"coherent" energy spread have well-established and precise definitions: "incoherent" means 
that the summation over different neutrino energies happens in the cross section of the total 
production - detection process whereas a summation over different neutrino energies in the 
amplitude is called a "coherent" energy spread.0 In the following we adhere to these defini- 
tions and refer the reader to section |VI| for a clarification of the notion of neutrino energy 
in the field-theoretical approach where the oscillating neutrinos occur in an inner line of the 
combined production - detection Feynman graph. Having different oscillation lengths in the 
process under discussion, then clearly the relevant coherence length is given by Jl7 



L™ h = minL$ fc . (1.4) 



Both kinds of neutrino energy spread, coherent and incoherent, lead to a loss of the oscillation 
pattern if L > L^ h and cannot be distinguished experimentally 0. Apart from the condition 
L < other conditions have to be fulfilled for the oscillation pattern to be present which 
do not depend on L ]3|, [13 , [14 . 



In this paper we will use the LSND experiment P| with the process 

/i + — > e + + v e + v e + p — > n + e + (1.5) 

as a model for our investigation for two reasons: first of all, the neutrino source (/x + ) is 
unstable and we want to extend the field-theoretical approach of Ref. [TB|] by taking into 
account the finite lifetime of the source; secondly, there is a claim made in Ref. |19| that in 



the LSND experiment the condition for coherence is not fulfilled. In the following, we will 
discuss in detail the effects of 



The notion of a neutrino wave packet is synonymous with the presence of a coherent neutrino 
energy spread. 

2 Hence, whether the summation over neutrino mass eigenstates or the summation over neutrino 
energies is concerned, "coherent" refers to a summation in the amplitude whereas "incoherent" 
refers to a summation over squares of amplitudes. 
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1. the quantum-mechanical uncertainties of momentum and energy of the initial particles 
involved in the production and detection processes represented by the widths of their 
respective wave packets or stationary states, 

2. the finite lifetime of the neutrino source particle, 

3. the uncertainties in the measurements of energies and momenta of the particles in the 
final state of neutrino production and, in particular, of the detection process. 

We are not able to take into account the interaction of the neutrino source particle (the fi + 
in our case) or the interaction of particles in the final state of the source process (in our case 
the positron originating from the fi + decay) with the matter background in which the source 
particle is generated (in LSND this background is water) in the field-theoretical treatment. 
We will only comment on the second of these effects in the last section of the paper. In the 
wave packet approach it is said that the interaction of the particles in the source process 
with the matter background interrupt neutrino emission and estimates of these effects are 
used to determine the "size of the neutrino wave packet" [Q. 

To include the finite lifetime of the neutrino source of the process ( p..5|) we combine 
field theory with the Weisskopf-Wigner approximation |20| in section O. In section [TIT] we 



calculate the amplitude for the reaction ( |1.5| ) by taking into account that the distance L 
between the source and the detector is macroscopic. This is achieved by using a theorem 



proved in Ref. |T3[ and an integral discussed in detail in the appendix of the present paper. 
In section |IV| we derive conditions for the existence of neutrino oscillations independent of 
L, while in section [V] we study some aspect of the cross section of the total production - 
detection process concerning the finite lifetime of the source. All the conditions for neutrino 
oscillations obtained in sections [TV] and[V] - whether dependent on L or not - are discussed in 



detail in section VI where we also study the problem raised in Ref. [IS] and try to elucidate 



the nature of the neutrino energy spread. 



II. PERTURBATION THEORY AND WEISSKOPF WIGNER 

APPROXIMATION 

To fix the notation we shortly repeat the basics of time- dependent perturbation theory. 
We consider a system described by the Hamiltonian H = H + Hi where H and Hi are not 
not explicitly time-dependent. The eigenvalues and eigenvectors of H will be denoted as in 
the relation Ho<f)j = Ej<pj where {<j)j} is a complete orthonormal system of states. For an 
arbitrary state if) = . Cj(0)(f)j at t = 0, the Schrodinger Equation gives the time evolution 
ip{t) = J2j c i(^)0j e_l£ ' : '* where the components Cj(t) obey the relations 

k k 

and 

H lM (t) = e iHot Hi(0)e- iHot (2.2) 
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defines the interaction Hamiltonian in the interaction picture. 

Let us now study the fi + decay and the "subsequent" detection of v e by u e + p — > e + + n. 
The interaction Hamiltonian is given by 

H 1 = H+ + H S +H+ + H D (2.3) 

where the indices S and D denote source and detection, respectively, and Hg and Hp are 
given by the Hamiltonian densities 

G 

T-Cs = T^ 7 ^ 1 ~ 75 ^ s ^Y(l - 7s)e , 

= -^=cos^cz/d7a(1 -7 5 )en7 A (l - fiU7s)p (2.4) 

describing muon and neutron decay, respecitively, and 

vs = U^Vj , v D = U ej Uj . (2.5) 

Actually, v e in Eq.( |2.4|) should be replaced by U*ji>j, however, this has no effect on the final 
result for neutrinos much lighter than the mass of the muon. The Hamiltonians with the 
superscript + are the Hermitian conjugates of those which carry the minus sign. 

Let us now sketch how to incorporate the finite muon lifetime in perturbation theory 
20| . To this end we take into account Hp, only when it occurs together with Cj(t) (see the 



initial conditions (|2.?1)) but take Hs in all instances. With this proviso the following states 
and coefficients are involved in perturbation theory: 

initial state p Cj, fa 

intermediate states e + ,i , e ,u s ; p <-> c'j, </>'■ 

e + ,n,u D <- 4, <\>" k (2.6) 

e + ,v e ,v s ; e + ,n,v D «-> c'^ k , ^ k 

final state e + ,z/ e ; e + ,n <-> c/, 0/ 

The initial conditions for the coefficients are given by 

q(0) = 1 and C/ (0)=c;.(0) = 4'(0) = c;^(0) = 0. (2.7) 

The differential equations for the coefficients are 

ic^-^'MMHsMMti), (2-8) 

3 

= E c K*)(^/I^Aint(O0i) + E^W^/l^ntW^') , (2-9) 

2c;.(t)^Q(t)m^+ nt (t)0 4 ), (2.io) 

i4(t) ~ Ci (t)(^|F+ int (t)0,,) + X: C*(*)(^l^int(*)^> , (2-11) 

^ ^(tX^I^Ct)^ . (2.12) 



With these approximations we get a closed system for q, c'j, c'l, c"'^ k . If we insert Eq. ( |2.12| ) 
into Eq.( p.llj ) we arrive at 



However, looking at the intermediate states ( |2.6|) , we see that the equations 
WL k \Ht M mi) = mHt M {V)<t>i) and E>> - E f » m = Et - E' 



(2.14) 



hold trivially because in the first matrix element e + , n, z/q and in the second second matrix 
element the proton are only spectators. Inserting Eq. ([2.14 ) into Eq. ( |2. 13|) then with a partial 
integration the second term of Eq. Q2.13Q is written as 



J,(Ei—E'.)(t—t') 



-m - e' 



k4V 







t i(Ei-E'Mt-t') 
dt' 6 „„ _,. cg(fl 



where 



[2.15) 



(2.16) 



We neglect now the term with c" in Eq. ( |2.15[ ) because it is of higher order and replace Ei — E 1 - 
by Ei — Ej + ie (e I 0) in Eq. fl2.15|) to have a well-defined expression. Since c'^(0) = we 
obtain 



3 

E A > P 



[i^-E'^t-t') 

-i(Ei-E!:) + e 



4(t' 



J 



^ 'Ei-E'j+ie 



Et-E' 



— Z7T 



- 4) <£(*) = (m - l -T)4(t) (2.17) 



where P denotes the Cauchy's principal value and T the total decay width of the muon. We 
neglect in the following the energy shift AEj or we can think it being already incorporated 
in the muon mass. Hence we get 



With similar arguments one obtains |2C| 

Ci (t) ~ e~5 r * 

and therefore 

Jo 

Using Eqs.(|2.9|) and ( |2.10| ) we get the final result 



(2.18) 



(2.19) 



(2.20) 
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c f (t) ~ (-if [ dti [ dt 2 
Jo Jo 

x (M (H+.J tl )Ht M (h)e~^ + H^Me-^H^ih)) 0,) . (2.21) 

This formula corresponds to the intuitive expectation. Apart from starting the time inte- 
gration at the inital time U — instead of £j = — oo we have the usual time-ordered product 
with the finite lifetime incorporated in the exponentials. 

III. THE AMPLITUDE 



With Eq. Q2.21p , the Hamiltonian densities ( p.4|) and Eq.( |2.2[) we can write for the am- 



plitude of the process ( |1.5| ) in the limit t — > oo 

A= (-i) 2 (Mp'u),4(p'es)i e D(PeD),n{p' n )\ 



xT 



dti J d 3 x, J dt 2 J d\ 2 n% hit {x 1 )e-^H + D ^{x 2 ) 



(3.1) 



where T is the time-ordering symbol. We assume that the muon fi + and the proton p are 
localized at the coordinates x s and x^, respectively. We imagine the proton being the 
nucleus of a hydrogen atom and bound in a molecule. Therefore we assume the proton 
state as stationary whereas the decaying muon will be described by a free wave packet with 
an average momentum equal to zero. This situation corresponds to the LSND experiment 
where the /i + is assumed to decay at rest. Since neutrino production and detection are 
localized at xs and xd, respectively, the spinors of the initial particles can be written in 
coordinate space as 

^ P {x)=^ p {x-x D )e- iE * t (3.2) 

and 

M*) = J MP) e-** 3 -**™ x e «**s , (3.3) 

respectively, with E^{p) = v/m^ + P 2 - The function ip p (y) is peaked at y = and the wave 
packet ipnip) in momentum space is peaked around the average momentum (p) = 0. The 
final particles will be described by plane waves. 

With the neutrino propagators of the mass eigenstate neutrinos 

WM^-W) = if Ifj^^^hr''^' < 3 - 4 > 

we obtain the amplitude 
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2 a]^a r _fp_ r M . f^, fV /■ 



* = y <^ x * y ^ y * y «** y (ar) 

x exp{i(p^+p; s ).xi + iO/ n +p^)-x 2 }exp{i(p.f 1 - E^t, - p-xs^e'^e'^ 
X^(P)7,(1 " 75)^ £ ^ ^^Iri^ ^Va - ls)ve{ P ' eD ) 

U 'it: T it 

XJs(pl,Pes) Un{p' n )l\{l ~ #A7s)Vv(^2 - ^d) (3.5) 

with 

WuM = ^M)Ya - Ileitis) ■ (3-6) 

We start with the integration over t x where we have to calculate the integral 

poo 1 

/ e - i{qo -El-E> eS+ E»mti e -ir tl dt = ^ ( g ?) 

Jo *(go-^-^ 5 + ^(p)) + ! r 

For the integration over t 2 we use the relation 



lim / e iEt2 dt 2 = iP[-)+ ir5(E) . (3.8) 



Hence a factor 



(3.9) 



appears in the amplitude. Furthermore, in the integration over x 2 we use the relation 

(2tt)- 3 / 2 y d 3 xe~ ils f(x + b) = e iU f(k) (3.10) 

where / is the Fourier transform of /. The integration over x\ is again trivial leading to the 
delta function 

(27r) 3 5(q-pl-p' eS +p). (3.11) 

Thus we obtain 

, G 2 F cos <d c . f ,s f d 4 q r ■/->! ->/ —\ —> i 

A= 1 J dp J j—^exp{-ip-x s -t{Pn+PeD + q)-x D } 



x(2n) 3 5(q-pl-p!,s + p)- 



i(q + E s ) + ±r 



i p(-±^-) + «6( go + E D ) 



x^(p)7p(i - 75)£f w - * o ' £C/r A (i - K)Ve(p' eD ) 

y '/^j "T" *t 

X-^GCPes) «n(K)7A(l - <M7s) Vv(<f + Pn + p'eo) (3.12) 

where we have defined 
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E s = E„(p) - E' v - E' eS and E D = E' n + E' eD - E p . (3.13) 

The integration over p can easily be carried out because of the delta function (|3.11|) and 
leads to the amplitude 

G 2 F cos -de 



A 



2^' 



-iq-L 



X 



i(q + E s ) + IT 



x^u(pi - g)7p(i - 75)^ 



7r5(g + ^D) 



1 



9o + 



« g2 



+ ie ' 1 



U* ejl \l- lh )v e (p' eD ) 



X JsiPu, Pes) U(p' n hx(l ~ 9Al^p{q + P2) 



where 



Pl=fil + PeS , P2=Pn+ p'eD and L 



x D -x S - 
E„ 



(3.14) 



(3.15) 



-q + pJ-El-E'e 



eSi 



Note that as a consequence of the integration over p we have Es — ^ 
i.e., is now a function of q. 

Now only the integration over q remains. Since we have a delta function of go within the 
brackets, the integration of the first of the two terms of the amplitude is trivial. We will 
show in the appendix that the go integration in the second term, which contains Cauchy's 
principal value, leads in the limit of a macroscopic distance L to the same result. In other 
words, in the limit of macroscopic L we have simply 27r<5(go + Ed) from the ti integration. 
In this limit we can apply a theorem proved in Ref. |13[ to perform the d 3 q integration and 
calculate the leading term of the amplitude for large L: 



A°° = J2 U ^ U e^ L A 



00 
3 



G 2 F cos$ c 2vr 2 



i{E Sj - e d ) + ±r 



x Vv(Pi + 1o l )7p(1 - 7b)(- hi + m j)7A(l - lb)v e {p' eD ) 

xJsiP'viP'esMp'J^i 1 ~ 9Al 5 )M-q/+P2) 



(3.16) 



where the definition of A°° is obvious, / is the unit vector pointing from the neutrino source 
to the detection point, kj are the momenta of the intermediate neutrinos and 



kj — 



E D 
qjl 



Qi 



Note that 



E S j = E^qjl+p l )-E' v -E' e 



eS 



(3.17) 



(3.18) 



because by virtue of the theorem in Ref. |T3| for each j the vector q has to be replaced by 
—qjl- The irrelevant phase factor occurring in the first line of Eq. ( |3.14)) has been dropped 
in A°°. 
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IV. COHERENCE CONDITIONS INDEPENDENT OF L 



Inspecting Eq. ( |3.16 ) it is evident that oscillations involving mh — m 2 k can only take place 
if 1,0 

\Qj ~ Qk\ < cr s and - q k \ < a D (4.1) 

where as and <td are the widths of ij)^ and ipp, respectively. We call conditions ( |4.1| ) ampli- 
tude coherence conditions (ACC). If they are not fulfilled, either by the source wave packet 
or the detector wave packet, then .4°° x A'j? ~ (j ^ k) which means that the term labelled 
by jk is suppressed in Eq. (|1 . 1|) . In the ultrarelativistic limit Eq. ( ^4.1|) is rephrased as 

Affl-L 

7/ £ °s,d ■ (4.2) 

Denoting by a x s,xD the widths of the wave functions Vv( x ) an d 4>p(. x ) m coordinate space, 
respectively, then, with the uncertainty relations a x s, x d &s,d > 1/2, Eq.(|4.2j) is rewritten as 

v*s,*n < (4.3) 

where we have made the identification E v = Ed (see Eqs.( [3.16|) and (|3.17| )). 
The amplitude Aj° contains the factor 

(4.4) 



i(E Sj - E D ) + |r 

which leads to a condition analogous to Eq.( |4.1| ) for neutrino oscillations to take place: 

\E Sj -E Sk \< l -T. (4.5) 

In the following this condition will be called source wave packet - finite lifetime condition 
(SFC). In the ultrarelativistic limit and with as <C we obtain 

/\ty) 2 _ 

E Sj - E Sk « ^(E D + l-pi) (4.6) 

and, assuming that ACC holds, we observe that \E D l + Pi\ < a s is valid (see the argument 
of ifjp in .4°° (|3.16|) ) and therefore Eq.( |4.5|) is rewritten as 



Am 2 jk a s 



m„E 



< r . (4.7) 



Defining At> M = as/ as a measure for the spread in velocity of the muon wave packet and 
with the muon lifetime r M = 1/T we can interpret Eq. (^4.7| ) as 

A W <i-L°r- (4.8) 
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V. THE COHERENCE LENGTH DUE TO THE FINITE fx LIFETIME 



Having performed all the integrations in the amplitude in the limit L — ► oo, we will 
discuss some aspect of the integration in the cross section. There we have integrations of 
the form d 3 p'/2E' for each particle in the final state, i.e., in our case u e and in the 
source process and n and in the detector process. In general these integrations cannot 
be performed without knowledge of the source and detector wave functions. However, for 

T < CT S ,D (5.1) 

the factors 

{{i{E Sj - E D ) + r/2) x {-i{E Sk - E D ) + T/2)}- 1 (5.2) 

in the cross section can be used to apply Cauchy's Theorem in order to obtain the coherence 
length associated with the finite muon lifetime. We assume that the ACC and SFC are valid 
and integrating over momenta of the final state of the detector leads to an integration in the 
variable Ed over a particular interval containing E$j Vj such that the length AEd of this 
interval fulfills T <C AEd <C <Js,d- This allows to write Ed = Ed + e where Ed denotes the 
central value of the interval which we define as 

E D = (E Sj ) (5.3) 

where (E$j) denotes the mean value of the E$j and the integration variable e varies over 
|e| < AEd/2. From the SFC it follows that Ed ~ Egj Vj which is exact up to terms of 
order T. With Am 2 - k > and qj — q% ~ —Am 2 - k /2E D we observe that the e integration 
over the interval on the real axis can be closed via a half-circle in the upper half-plane since 
the radius of the half-circle is much larger than V and therefore the factors ( |5.2|) make the 
part of the path in the upper half-plane negligible in the integral. Because of Eq.( [5.1|) and 
AEd ^ &s,d this integration does not affect the source and detector wave functions to a 
good approximation. Then Cauchy's Theorem states that the result of the integration is 
given by making the replacement 

E D - E^ = E D + l -T (5.4) 

(0) 



in the absolute square of ^4°° ( |3.16|) . Inserting E D into the exponential —iAm? k L/2ED we 



see that the cross section contains the damping factor 

Am%T 



exp(-^^L) . (0.0) 



As we will discuss in the next section the requirement ( |5.1|) is very likely to be fulfilled for 
LSND with the decay width of the muon being r ~ 3 x 10 -16 MeV. Ed can still be thought 
of as being identical with Ed (|3.13|) for all practical purposes because this T is so small that 
AEd can be chosen smaller than any achievable accuracy for Ed in a real experiment. 
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VI. DISCUSSION 



The characteristics of the field- theoretical approach: In this paper we have used the 
field-theoretical approach to discuss neutrino oscillations as we have done in Ref. ||13|| . In 
this approach the whole process of neutrino production and detection is represented by a 
single Feynman graph such that the oscillating neutrinos are associated with the inner line of 
the graph. However, due to the macroscopic distance between neutrino source and detector 
this inner line is on-shell (see Eq. fl3.16| )) for each neutrino with definite mass according to 
a theorem proven in Ref. |13| . In the present paper we have incorporated the finite lifetime 
of the neutrino source which is given in our concrete example of the LSND experiment (see 
Eq.( |1.5|) ) by a positively charged muon whose decay is responsible for the neutrinos with 
which the experiment is performed. The finite lifetime of the neutrino source prevented us 
from using ordinary perturbation theory with an initial time ti = — oo. Instead we took 
advantage of the Weisskopf-Wigner approximation which allowed us to take ti — and to 
combine in this way the finite muon lifetime with perturbation theory. 

With the help of the above-mentioned theorem all integrations in the amplitude of the 
combined production - detection process could be performed in the asymptotic limit L — > oo. 
From the requirement Af ~ A™ ( p,16|) we derived the amplitude coherence conditions 
(ACC) (|Op, O and the source wave packet - finite liftime condition (SFC) (fO|), (gTTJ). 
For a given mass-squared difference Am 2 fc these three conditions are the prerequisites for 
neutrino oscillations. If they are not fulfilled, the term with exp(—iAm^ k L/2E u ) is sup- 
pressed in the oscillation probability which means that no neutrino oscillations with respect 
to Am 2 fc are possible. Here we have identified the neutrino energy E v with E D ( EU3D which 
is justified in view of the definition of qj Q3.17 ) occurring in the exponentials e iqjL of A°° 
( |3.16| ). ACC and SFC are both independent of L, therefore no coherence lengths are associ- 
ated with them. In coordinate space Eq.( |4.1| ) simply means that the oscillation length must 
be larger than the widths of the production and detection wave functions (see Eq.( [4.3|) ). The 
ACC were among the main results of Refs. || and |13[ and an analogous condition in the 



framework of the wave packet approach was recently emphasized in Ref. [12]. The SFC (|4.5| 



says that the spreading of the muon wave function?] during the muon lifetime should be less 
than the oscillation length in order not to wash out neutrino oscillations (see Eq.( |4.8| )). In 
Ref. |TB[] this condition was not found because it was assumed that the source wave function 
is stationary. Clearly, in such a case the energy of the neutrino source is fixed and does not 
depend on rrij and, consequently, no SFC is present. 

Discussion of the LSND experiment: Coming back to the LSND experiment, as rep- 
resents the momentum spread of the stopped muon and an estimate of it is given by 
a s ^ 0.01 MeV |2T|. For the detector proton bound in CH2 groups (mineral oil) [18| 



it is reasonable to assume that in coordinate space its wave function has a spread of 
around 1 A and consequently o\d ~ 2 x 10~ 3 MeV. Dropping now the indices of the mass- 
squared difference, with representative values Am 2 ~ 1 eV 2 and E u ~ 30 MeV we obtain 
Am 2 /2E U ~ 10~ 14 MeV and we conclude that the amplitude coherence conditions are very 



The muon wave function is non-stationary. 
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well fulfilled in the LSND experiment. Performing an analogous estimate with the SFC we 
get Am 2 a s /m ll E u < 3 x 1CT 18 MeV < T ps 3 x 1(T 16 MeV with the numbers we used 
before for the ACC. From this result we conclude that also the SFC is valid in the context 
of the LSND experiment, though, surprisingly, the margin is only given by two orders of 
magnitude. Note that for the ACC this margin is eleven orders of magnitude. 

In Ref. |19] it was found that Am 2 /E u < T should hold for coherent neutrino oscillations. 
This condition is not fulfilled for LSND. In this paper with the field-theoretical treatment 
we could not recover this condition which looks like the first ACC (|4.2|) with as replaced 



by T. In this context we want to stress the following: in the process ( p..5| ) there are three 
different lengths, namely the sizes of the wave packets of the neutrino source of order 1 /as 
and of the detector particle of order l/a^, where as and a^ are the widths of ip^ and ip p , 
respectively, and the length 1/r associated with the finite lifetime of the source. Each of 
these lengths is uniquely defined and the function of each of them uniquely determined in the 
field-theoretical framework. However, in the wave packet approach the distinction between 
as, &d and V is not so clear and each of the three lengths could be associated with the size 
of the neutrino wave packet and possibly lead to erroneous conclusions. 

The characteristics of neutrino oscillations: We want to emphasize that in the field- 
theoretical approach the notion of a neutrino wave packet does not exist and the questions 
whether the neutrino energy or neutrino momentum is fixed or both can vary do only make 
sense in connection with the processes of production and detection. This is because only 
parameters associated with particles of the exterior legs of the Feyman graph, i.e, with those 
particles which are manipulated in the experiment, determine the neutrino oscillations. Let 
us notice that, for fixed momenta of the final state particles of the production and the 
detection processes, with Eq.( |3.16| ) the oscillation probability has the form 



(6.1) 



and we can imagine that the phase factors e iqjL represent the plane waves of the different 
neutrinos mass eigenstates. We will use this fact to compare the wave packet approach with 
the result of this paper. We arrive at the following characteristics of neutrino oscillations in 
our field-theoretical approach: 

1. We have chosen the detector wave function, i.e., the wave function of the proton, to 
be a bound state and therefore the detector wave function does not spread in time. 
This looks physically very reasonable to us meaning that the detector is always on and 
waiting (see, however, Ref. for a discussion of source and detector with a temporal 



resolution). As a consequence we have E u = E D (see Eqs.( |3.13"D , ( |3.16|) and ( |3.17| )) 



and Eq.(|6.1[) suggests that neutrino oscillations take place between neutrinos with the 



same energy but with different momenta qj ||13 



The identification E v = Ed allows to determine the neutrino energy with arbitrary 
accuracy by measuring the energies of the particles in the final state of the detector 
process, in our case the neutron and the detector positron, with arbitrary accuracy. 
Therefore, one can limit the averaging over the neutrino energy to an arbitrarily small 
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interval - of course, in practice at the expense of the number of events - and make 
any coherence length arbitrarily long by performing only detector manipulations (see 
Eq.( |1.3| )). This so-called restauration of coherence is trivial in our approach. It agrees 
with observations in the wave packet treatment [7,|11]. We find no upper limit to the 



coherence lengths in contrast to Ref. [0 which is due to the fact that we assume 
that the detector is sensitive to energies and momenta of the particles produced in the 
detector process whereas in Ref. [I2| it is assumed that the detector measures neutrino 
wave packets of a certain width. 

3. Taking the field-theoretical approach to neutrino oscillations seriously, assuming that 
the detector particle is initially in a stationary state and that the observation of par- 
ticles associated with the neutrino detection is done by energy and momentum mea- 
surements, we come to inevitable conclusion that there are no neutrino wave packets 
in neutrino oscillations because all summations over neutrino energies happen in the 
cross section and are thus incoherent summations.^ Our assumptions include the cases 
that some particles are not detected at all or that cuts in energy and momentum are 
made. In addition, any further measurements of observables commuting with the en- 
ergy and momentum operators performed by the detector do not change our conclusion 
[[?]]. We think that our conclusion is correct for realistic experiments. If one assumes 
instead that, with respect to the particles in the final state of the detection process, 
the detector is sensitive to wave packets of some form then our conclusion is not valid. 
However, we do not know if such a detector exists. 

The coherence lengths: Let us now assume that the ACC and SFC hold and study 
the effects of different energy and/or momentum averaging mechanisms which all lead to 
specific coherence lengths. In the light of the above discussion all coherence lengths result 
from incoherent neutrino energy spreads. There are three types of coherence lengths |7j 
associated with our neutrino production and detection processes (|1.5|) : 

A. L c ^ h due to the finite muon lifetime, 

B. L c g h due to the interruption of the neutrino emission because of collisions of the source 
positron with the background and 

C. L£? h due to the neutrino energy spread introduced by the usual imperfect energy 
measurements of the particles observed with a realistic detector. 

According to Eq. (|5.5| ) the first of the coherence lengths is given by 



Though in this paper we consider a neutrino source at rest, this conclusion and points 1, 2 are 
also valid for accelerator neutrinos because the arguments leading to it depend on the neutrino 
detection process and not on the production. Formulated in another way, it is the detector which 
through its properties determines the nature of the neutrino energy spread. 
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This coherence length appears if the energy averaging amounts to a summation of Ed over 
an interval much larger than T (see section With 1/T corresponding to 659 m and 
Am 2 ~ 1 eV 2 , E v ~ 30 MeV this coherence length is around 200 light years and completely 
irrelevant for the LSND experiment^ and, therefore, apart from its effect in the SFC, the 
finite lifetime of the muon could have been neglected as was done in Ref. with the 
lifetime of the neutrino source nucleus. However, the main point in the investigation of L™ h 
was rather to see how it emerges in the field-theoretical treatment. It is interesting to note 
that L c ^ h enters into the cross section through exp(— L/L c ^ h ). We have obtained this form 
of the damping factor for T <C 05,00. If this condition is not fulfilled the damping factor 
depends on the form of the functions ip^, ip p . In the approach using Gaussian wave packets a 
corresponding damping factor has the form exp{ — (L/ L coh ) 2 } |12|, however, it is not obvious 
how to compare this factor with the previous one and which coherence length is meant with 
L coh . 

L°g h is not included in our treatment because we do not know how to deal with random 
collisions of the source positron with the matter background and the associated heuristic co- 
herence length ( |1.3| ) in the field-theoretical approach. So here we only repeat the arguments 
found in the literature about the neutrino wave packet approach. There one estimates the 
mean free path of the positron from the fi + decay in the matter background where fi + is 



produced. In the LSND experiment this background consists of water |18| and according 
to the rule of thumb presented in Ref. || this mean free path of the positron is of the or- 
der of centimeters. Then in the heuristic approach one would estimate the size £b of the 
neutrino wave packet to be of the same order of magnitude @. Thus, adopting the wave 
packet approach, one would estimate L c ^ h ~ E 2 £b/ 'Am 2 which is something like four orders 
of magnitude smaller than L™ h , but still astronomical, making the previous consideration 
of the coherence length originating from the finite muon lifetime (or a wave packet size of 
order £a ~ cr^ = 659 m) spurious. 

The coherence length |7J Lg> h = L osc E v /AE Ul where AE V comes from the inability to 
measure E v = Ed better than at a certain realistic experimental accuracy, is the only relevant 
coherence length in practice |[L7|| . In the LSND experiment [jl8[] the detector positron and 
neutron are detected in coincidence and, in addition, the energy of the positron is measured. 
We guess that AE V is of the order of 5 MeV [18], thus Lg? h is probably not more than several 
times the oscillation length. This is many orders of magnitude smaller than the astronomical 
coherence lengths L c ^\. 

Summary: Let us summarize the main points. We have confined ourselves to situations 
where the neutrino source is at rest, thus the present investigation is not straightforwardly 
applicable to high energy neutrinos (see, however, Ref. WQ\ for a field-theoretical discussion 
of accelerator neutrino experiments). We have assumed that the wave function of the de- 
tector particle is stationary. Then the field-theoretical approach to neutrino oscillations is 
completely static and there are no explicit time averages necessary as in the wave packet 



5 Note that even this coherence length could theoretically be overcome by measuring Ed with a 
precision better than V and chosing events with Ed in a given interval of size AEd smaller than 
r (compare section [V] and point 2 in this section). 
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approach. With the field-theoretical method we have clarified, using the model reaction 
(1.5), the roles of the widths as, <Jd and T in neutrino oscillations where these widths are 



associated with the source, the detector and the finite lifetime of the source, respectively. 
To check the validity of the ACC (|4.2|) and the SFC ( f4.7|) in a real experiment concrete 
values of as and a® have to be chosen. Making a plausible guess for an and using as < 0.01 
MeV |2]| in the case of LSND, the ACC are very well fulfilled and also the SFC seems to 
hold safely. Finally, if our method is a correct approach to neutrino oscillations then, in 
experiments with the realistic detector properties assumed in this paper, oscillations take 
place between neutrino mass eigenstates with the same energy but different momenta, there 
are no neutrino wave packets and the coherence length in neutrino oscillations results from 
an incoherent neutrino energy spread. 
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APPENDIX A: INTEGRATION OVER q 



We consider the integrals 



<'</" i—TT- P ( — T~ ) .2 J'^L I T. ( Al ) 



i(q + E s ) + ±T \q + E D J q 2 - q 2 - mj + ie 

—oo 

where n = 0, 1 and P denotes Cauchy's principal value. We will calculate this integral with 
the help of the residue calculus which gives the formula 

dxP (_L\f( x ) = Uf + M dx _^ /(x ) (A2) 
\x-x J 2 \J Cr JcJ x-x 

for a function / which is analytic along the real axis. The paths Cr and Cl lead along the 
real axis except close to x = Xq where the point x = Xq is circumvented to its right or to its 
left in the complex plane, respectively. In our case Xo = —Ed and / has three poles of first 
order at 



qf = -E s + l -T , q^ = ^q 2 + m) - it and q^ = - ^q 2 + m 2 + ie . (A3) 
Since we have only one pole below the real axis we close the contour below. Then we obtain 

£S , . 



In the second part of Eq.(|A|) we use that -E^ > and perform the limit e — > without 
getting a singular integrand. The integral I n depends on q and appears in the amplitude in 
the following way: 



A = J d 3 q<5>(q)I(q)e- ic ? L (A5) 
where $ can be read off from Eq.( |3.14|) . For the second term of I n we can use Lemma 3 



of Ref. [13] to show that it decreases like L~ 2 for L —>■ oo. Since we are only interested 
in the leading term oc L~ x of the amplitude for large L we can neglect the contribution of 
the second term of I n . It is then easy to show that the first term in the integral (|A4|) gives 



exactly the contribution to the amplitude ^4°° ( ^.16| ) as the term with ir5(qo + Ed). 
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